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Electronic state and the Verwey transition in magnetite (Fe304) are studied using a spinless 
three-band Hubbard model for 3d electrons on the B sites with the Hartree-Fock approximation 
and the exact diagonalisation method. Complex-orbital, e.g., -^[|za;) + i|j/z)], ordered (COO) 
states having noncollinear orbital moments ~ 0.4 [ib on the B sites are obtained with the cubic 
lattice structure of the high-temperature phase. Fhe COO state is a novel form of magnetic 
ordering within the orbital degree of freedom. It arises from the formation of Hund's second 
rule states of spinless pseudo-d molecular orbitals in the Fe4 tetrahedral units of the B sites 
and ferromagnetic alignment of their fictitious orbital moments. A COO state with longer 
periodicity is obtained with pseudo-orthorhombic Pmca and Pmc2\ structures for the low- 
temperature phase. The state spontaneously lowers the crystal symmetry to the monoclinic and 
explains experimentally observed rhombohedral cell deformation and Jahn- Teller like distortion. 
From these findings, we consider that at the Verwey transition temperature, the COO state 
remaining to be short-range order impeded by dynamical lattice distortion in high temperature 
is developed into that with long-range order coupled with the monoclinic lattice distortion. 

KEYWORDS: magnetite, Fea04, the Verwey transition, complex-orbital order, charge order, noncollinear 
orbital moment, pyrochlore lattice, multiferroic 



1. Introduction 

Magnetite (FeaO^) is one of the most familiar miner- 
als for human beings. Known as a loadstone even in the 
Greek era, it has been attracting us for its magnetic prop- 
erties. In 1939, Verwey discovered an abrupt increase of 
resistivity in magnetite by two orders of magnitude on 
cooling it below Ty ~ 120 K. 1 Despite numerous exper- 
imental and theoretical studies, the mechanism of this 
Verwey transition has been controversial for more than 
half a century. 

Magnetite is a kind of spinel ferrites. It crystallises in 
the cubic inverse spinel structure at room temperature 
(see Fig. 1). One-third of Fe ions are on tetrahedral A 
sites with a Fe 3+ formal valence constituting a diamond 
lattice. The remaining two-thirds of Fe ions are on oc- 
tahedral B sites in a Fe 2 5+ formal average valence state 
forming a pyrochlore lattice. Below Tn ~ 860 K mag- 
netite is a ferrimagnet, where the magnetic moments of 
Fe ions on the A and B sites are antiparallel and the 
net magnetisation is ~ 4 /is per unit formula. The tetra- 
hedral (octahedral) oxygen coordination of the A (B) 
site Fe ion causes ligand field splitting of the 3d level 
into three-fold t2 g and two-fold e g levels. The 3c? orbitals 
with minority (majority) spin on the A (B) site are fully 
occupied and half of Fe ions on the B sites accommodate 
an extra electron with minority spin in their t2 g orbital 
in each. 

The Verwey transition is a first-order phase transi- 
tion accompanied by a lattice deformation from the cubic 
(Fd3m) symmetry. Below Ty, the cubic cell is elongated 
to nearly rhombohedral along the [Hl] c direction (the 
subscript 'c' denotes the cubic cell). The lattice sym- 
metry of the low-temperature phase is still unresolved 
problem. X-ray and neutron diffraction studies have 




Fig. 1. Illustration of crystal structure of Fe304 in the high- 
temperature cubic cell. The large open and hatched circles repre- 
sent Fe ions on the B and A sites, respectively and small circles 
denote oxygen ions. 



shown that the low-temperature structure has a mono- 
clinic Cc symmetry with a \f2a c x V2a c x 2a c crystallo- 
graphic supercell. Small triclinic distortion from the Cc 
structure has been inferred from the observation of the 
magnetoelectric effect 5,6 and a study of twins with the 
synchrotron radiation x-ray topography. 7 However, no 
evidence has been found for the breaking of the inversion 
symmetry in an infrared spectroscopic investigation. 8 

Verwey and co-workers have interpreted the transi- 
tion as an order-to-disorder transition of the Fe 2+ and 
Fe 3+ ions on the B sites; Fe 2+ and Fe 3+ layers stack al- 
ternately along the [001] c direction so as to minimise 
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Coulomb energy in their charge-ordering model. 9 Al- 
though intensive experimental and theoretical efforts 
have been made, the existence of the charge ordering in 
the low-temperature phase as well as its electronic states 
above and below Ty still remain matters of debate. 10-12 

Several theoretical models have been proposed to ex- 
plain the mechanism of the Verwey transition. Cullcn 
and Callen have discussed the transition using a spinlcss 
one-band Hubbard model within the Hartree approxi- 
mation. 13 They assumed that the trigonal crystal-field 
is strong enough to lift all the degeneracy of the t 2g or- 
bitals and only the low-lying a\ g orbital of every Fe ion 
on the B sites is taken into account. They found that 
a charge-ordered state is stabilised when Coulomb inter- 
action between the neighbouring sites exceeds a critical 
value. Mishra et al. have extended the one-band model 
to a three-band model considering the triple degeneracy 
of the t 2g orbitals. 14 Their results are similar to those 
of Cullen and Callen's except for the ferro-orbital order- 
ing, where only one kind of the t 2g orbitals (xy, yz or zx 
orbitals) is taken for all the occupied sites. 15 Based on 
the three-band model, Seo et al. have proposed a mech- 
anism of the Verwey transition with a bond dimerisation 
caused by the Peierls instability, where charge ordering 
is irrelevant to the transition. 15 Recently, the electronic 
state in the low-temperature phase has been investigated 
using local density approximation + Hubbard U (LDA 
+ V) calculations 16-19 based on the monoclinic lattice 
structure shown by Wright et al. 4 In the LDA + U cal- 
culations, a state with simultaneous charge and antiferro- 
orbital ordering is obtained. 

Large orbital moments at the B sites above and be- 
low Ty have been observed in recent measurements of Fe 
^2,3-edge x-ray magnetic circular dichroism (XMCD). 20 
However, more recent XMCD experiments have resulted 
in observing no net orbital moments. 21 On the other 
hand, orbital moment ~ 0.5 /xb has been also inferred 
from the magnetic Compton scattering experiments. 22 

In previous works, we have investigated a spinlcss 
three-band Hubbard model with the Hartree-Fock (HF) 
approximation and exact diagonalisation method and 
found complex-orbital ordered (COO) states with non- 
collinear orbital moments as the ground states within a 
realistic parameter range. 23,24 Here, the word 'complex'- 
orbital represents an orbital which is described as a linear 
combinations of the t 2g basis functions (xy, yz and zx 
orbitals) with complex number coefficients and has non- 
zero orbital moment. 

The purpose of this paper is to describe the forma- 
tion mechanism of the COO state and a possible COO 
state realised in the low-temperature phase and to dis- 
cuss the mechanism of the Verwey transition. We found 
that the COO state is a kind of magnetic state within 
the orbital degree of freedom. It does not originate from 
the mechanisms in the conventional orbital order, which 
can be described with the Kugel-Khomskii type model 25 
or in usual magnetic materials, where orbital moments 
are induced by ordered spin via the spin-orbit interac- 
tion. A COO state with longer periodicity is expected 
to be stabilised in the low-temperature phase owning to 
the strong coupling between the orbital polarisation and 



lattice distortion. The state explains experimentally ob- 
served rhombohedral cell deformation and Jahn- Teller 
like distortion. 

This paper is organised as follows. The details of the 
model Hamiltonian are presented in § 2. Charge and or- 
bital ordered states found in the HF calculations with 
the cubic structure (T > Ty) are described in § 3. The 
formation mechanism of the COO state is discussed ana- 
lytically and the stability of the state against the electron 
correlation effects is examined with the exact diagonal- 
isation method for a 16-site model in § 4. After short 
introduction on the lattice structure and related proper- 
ties in T < Ty, discussion on the COO state in connec- 
tion with lattice distortion and description on a COO 
state obtained from the HF calculations with pseudo- 
orthorhombic structure models for T < Ty arc provided 
in § 5. In § 6, electronic states above and below Ty and 
the mechanism of the Verwey transition are discussed. 
Finally, a summary is given in § 7. 

2. Three-Band Hubbard Model 

To discuss possible charge and orbital orderings in 
magnetite, we studied a spinless three-band Hubbard 
model for 3d electrons on the pyrochlore lattice of the 
B sites. The pyrochlore lattice is a network of corner- 
shared tetrahedra consisting of linear chains with six dif- 
ferent directions [1, ±1, 0] c , [1, 0, ±l] c and [0, 1, ±l] c (see 
Fig. 1). This model was originally introduced by Mishra 
et al. 14 as an extension of the one-band Hubbard model 
of Cullen and Callen. 13 In the three-band model, only t 2g 
electrons with minority spins at the B sites are considered 
on the assumption that the spin moments at the B sites 
are collinear, and the 3d orbitals with majority spins are 
fully occupied. We furthermore introduce effects of the 
spin-orbit interaction, a trigonal crystal-field and elec- 
tron hopping via adjacent O 2p orbitals. As shown later, 
electron hopping among the B sites via the neighbour- 
ing oxygen ions is indispensable for the complex-orbital 
orderings. 

The Hamiltonian for t 2g electrons with minority spins 
on the B sites is written as 

H = E ^(i^ + sW 
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Here t%j denotes hopping integral between orbital \x 
(~xy, yz or zx) at site i and orbital v at site j. c\ , 
and rii^ represent the creation, annihilation and num- 
ber operators of an electron on orbital fi at site i, re- 
spectively, (ij) stands for summation over the nearest- 
neighbour sites. U and V are the on-site and inter-site 
Coulomb energies, respectively. The first term in the 
square brackets on the third line describes the trigonal 
crystal field, where I is an I = 1 pseudo-angular mo- 
mentum operator, and the unit vector on is along the 
three- fold axis of site i. The trigonal field splits the t 2g 
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levels at a B site into a single a\ g and two-fold e g levels, 
and the level splitting is denoted as D(= e(e^) — e(a lg )) . 
The second term in the square brackets is the spin-orbit 
interaction with a coupling constant £ = 0.052 eV; the 
unit vector as represents the direction of the minority 
spins at the B sites. 

Hopping integrals between nearest-neighbour sites i 
and j can be describe with the Slater-Koster parameters 
tdda, tddir an d t dd s as 26 

= 3l 2 m 2 t dda + (I 2 + m 2 - 4l 2 m 2 )t ddv , 
t yz,zx _ 2,l mn 2 t dda + lm(l - 4:n 2 )t dd7r - sgn(/m) t pd , 

(2) 

where I, m and n denote direction cosines (I, m, n) along 
site i to site j direction and terms with t dd $ arc omitted. 
The additional term with t pd in t yz,zx describes electron 
hopping through the adjacent O 2p z orbital. There is the 
relation t^ v = t^, and one can obtain all other terms by 
interchanging the x, y and z coordinates. The relations 
tdda = — 63.2r~ 5 eV and — 34. lr~ 5 eV are assumed 
after Harrison, 26 where r (in units of A) represents the 
distance between sites i and j. 

The values of hopping integrals for the B 
sites along the [1,±1,0] C chains in the cubic 
structure are t*? ,xy = 3t ddtT /4 = -0.206 eV, 

t y Z ,y Z = t z*M = tM j 2 = 0Q74 cV 

tf z ' zx =t zx ' yz = ±(t dd j2-t pd ) = ±(0.074 cV-t pd ), and 
all other terms are zero. t pd corresponds to the hopping 
process: Fe 3+ (3d 5 ); 2 "(2p 6 ); Fe 2+ (3d 6 ) -> 
Fe 2+ (3(i 6 );0"(2 2 3 5 );Fc 2+ (3d 6 ) -» 
Fc 2+ (3d 6 ); 2 ~(2p 6 );Fe 3+ (3d 5 ) and the value is roughly 
estimated to be t pd = t 2 d7T /A ~ 0.36 eV by means of the 
second-order perturbation theory. Here, t pd7V ~ 0.6 eV 
is the Slater-Koster parameter between Fe 3d and O 2p 
orbitals, and A ~ 1 cV is the charge-transfer energy 
at the Fc 3+ ion site defined as the averaged energy 
difference of the 3d 6 with a hole on the ligand and the 
3d 5 configurations. 27 The value of t pd is comparable 
to the direct d-d hopping terms and not negligible. 
U = 4.0 eV is adopted. 28 We do not know the accurate 
value of D, and here, D = 0.25 eV is assumed; the D 
dependence of electronic state will be discussed in § 4.2. 

3. Electronic State in the Cubic Structure 

Electronic structure calculations for the high- 
temperature phase were performed within the HF ap- 
proximation assuming the cubic fee unit cell with four in- 
dependent B sites (see Fig. 1). The wave vectors through- 
out the first Brillouin zone of the fee lattice were sampled 
using the tetrahedron method. 29 In the HF approxima- 
tion, we have 

+ {c\^ jv ){c) vCill ). (3) 

Only the Hartree terms [the first three terms in eq. (3)] 
were considered in the previous studies; 14 ' 15 however, 
inclusion of the Fock terms [the last three terms in 
cq. (3)] are essential to keep the spherical symmetry of 



the Coulomb interaction. While orbital ordcrings with 
variety of orbitals described by linear combinations of the 
three t 2g orbitals are obtained with the Fock terms, lim- 
ited orderings with pure yz, zx or xy orbitals are found 
without the term. 

In Fig. 2, phase diagrams obtained with D = 0.25 eV 
assuming the directions of minority spins along direc- 
tions [lll] c (a), [110] c (b) and [001] c (c) are shown. For 
comparison, that calculated without the spin-orbit in- 
teraction (£ = 0) is also depicted in (d). In each pan- 
els, six different phases COO-I, COO-II, ROO-I, ROO- 
II, CD-ROO-I and CD-ROO-II can be seen, except for 
the additional small area of COO-III phase in Fig. 1(c). 
These phases can be classified into two groups whether 
ordered orbitals are described as linear combinations of 
the xy, yz and zx orbitals with real number coefficients 
(labelled with 'ROO') or complex number coefficients (la- 
belled with 'COO'). The CD-ROO-I and CD-ROO-II 
states are real-number orbital-ordered (ROO) states ac- 
companied by charge ordering where the layers of Fe 2+ 
and Fe 3+ ions are alternately piled up along the [001] c 
direction. For small values of t pd the ROO states are sta- 
ble, whereas the complex-number orbital-ordered (COO) 
states are favourable for large t pd values. The ROO states 
with charge ordering prevail with large values of V. Dif- 
ferences among the phase diagrams are rather limited; 
large variations can be found only in the boundary be- 
tween the COO-I and COO-II phases and the appearance 
of the COO-III phase with a s //[001] c . 

The dig orbital of every B site is occupied by 0.5 elec- 
tron in the ROO-I state, a similar state but with an 
additional ferro-orbital component (one of the xy, yz 
and zx orbitals) is stabilised in the ROO-II state. For 
D = eV, we found the ferro-orbital ordered state cor- 
responding to that discussed in the Hartree calculations 
in ref. 15 with small value of V, in the range of 0.05 eV 
< t pd < 0.15 eV. The ROO-II emerges as a result of 
mixture of the ferro-orbital ordering and the a\ g orbital 
ordered state that are strengthen by the trigonal fields 
(D = 0.25 eV). The CD-ROO-I and CD-ROO-II states 
are accompanied by the Verwey charge orderings, where 
the layers of Fe 2+ and Fe 3+ ions are alternately piled up 
along the [001] c direction. The a\ g orbital of every Fe 2+ 
ion site is occupied in the CD-ROO-I state. In the CD- 
ROO-II state, ^[Ij/ 2 ) + \ zx )\ orbitals are occupied in 
Fe 2 + ions on [110] c chains as the lobes of the orbitals are 
parallel to the Fe 2+ chains. The ROO-I and CD-ROO-I 
states found in the small t pd regime are essentially the 
same as those discussed by Cullen and Callen, 13 where 
the orbital degeneracy is lifted. 

Figure 3 shows the directions of the orbital moments 
in the COO-I, COO-II and COO-III states. The chiefly 
occupied t 2g orbitals (~ 0.5 electron per site) and their 
orbital moments are listed in Table I. The total orbital 
moments are along the [lll] c , [H0] c and [001] c direc- 
tions in the COO-I, COO-II and COO-III states, respec- 
tively. These COO states have the same periodicity to 
the crystal, but the orbital ordering spontaneously lower 
the space group symmetry from cubic (Fd3m) to rhom- 
bohcdral (R3m), orthorhombic (Imm2) and tetragonal 
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Fig. 2. Phase diagrams of the ground states for the high-temperature phase in V-t p d plane with D = 0.25 eV. In the calculations, the 
directions of minority spins are assumed along directions [lll]c (a), [110] c (b) and [001] c (c). The same but without the spin-orbit 
interaction is also shown in (d) for comparison. 'ROO' and 'COO' represent real- and complex- orbital ordered states, respectively, 
and 'CD' denotes the state with the Verwey charge ordering. Weak charge disproportionation occurs in the COO-I and COO-II states 
(for details see text), and the charge distributions in the ROO-I, ROO-II and COO-III states are uniform. 
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Fig. 3. Alignment of orbital moments in Fe404 cube for the COO-I (a), COO-II (b) and COO-III (c) states. The large and small 
circles signify Fe ions at the B sites and oxygen ions, respectively. The arrows on the four independent B sites labelled 1-4 denote 
the directions of the orbital moments. The grey (black) arrows represent electron-poor (electron-rich) sites in (b). The magnetic easy 
axes of the COO-I, COO-II and COO-III states are oriented along the [lll] c , [110] c and [001] c directions, respectively. Details of the 
electronic states on each of the COO states are listed in Table I. 



( /42m ) , respectively. Because of noncollinear alignment 
of the moments, the size of total orbital moments are 
rather small. For instance, in the COO-1 state, each B 
site has orbital moment of about 0.5 /ib- However, the 
total orbital moment is only ~ 0.20 /iB per B site, since 
the orbital moment at site 1 is in the opposite direc- 
tion to the sum of those at sites 2, 3 and 4, and this 
almost cancel out the net orbital moment of each Fe 4 4 
cube. As was mentioned before, the spin-orbit interac- 
tion is not essential ingredient for the stability of the 



COO states. This becomes evident when the results ob- 
tained with [see Figs. 2(a)-(c)] and without [see Fig. 2(d)] 
the spin-orbit interaction are compared, where the phase 
boundaries between the COO and ROO states are almost 
unchanged among the diagrams. On the other hand, rel- 
ative stability among the three COO states are affected 
by the direction of the spin moment, since the differences 
in the total energies among the three are only ~0.01 eV 
per site. Indeed, as clearly seen in Figs. 2(a)- (d), each 
of the COO states is more stabilised when ag is along 
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Tabic I. Details of the COO-I, COO-II and COO-III states obtained with the high-temperature cubic phase. The direction of trigonal 
axis ctj, dominantly occupied orbital (with ~ 0.5 electron) and its orbital moment (in units of fig) on each of four independent B sites 
(see Fig. 3) are listed for each COO state. 
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the direction of its total orbital moment as compared the 
other directions. 

The charge disproportionation of the COO states is 
rather moderate as compared to the CD-ROO-I or CD- 
ROO-II states. The charge disproportionation of the 
COO-II state approaches uniform with decreasing value 
of V and increasing t p d, and the differences in the oc- 
cupation number between electron-poor and electron- 
rich sites are from 0.05 to 0.29 in the phase diagram 
of Fig. 2(b). On the other hand, the COO-I state is con- 
siderably uniform, and the differences in the occupation 
number among sites are at most 0.09 in the phase dia- 
gram of Fig. 2(a). 

The band structure and the density of states (DOS) for 
the COO-I state are shown in Fig. 4. In our calculations, 
the wave vectors throughout the first Brillouin zone of 
the fee lattice were sampled. Note that orbital orderings 
other than the ROO-I state spontaneously break the cu- 
bic symmetry, and thus some symmetry points in the 
Brillouin zone of the fee lattice become inequivalcnt. For 
example, the L point at (|, |, 5) and that at (-5, 5, 5) 
are not equivalent in the COO-I state whose trigonal axis 
is along [11 l] c shown in Fig. 3(a). Within the mean-field 
approximation, only a particular orbital can be occupied 
on each sites to reduce Coulomb repulsion energy, in the 
large U limit. The four lowest band in the figure, in- 
deed, correspond to the states mainly consisting from 
the orbitals listed in Table I. The low symmetries in the 
occupied orbitals and the large components of hopping 
matrices which hybridise different kinds of orbitals, e.g., 
tfj' zx make these bands repulsive in each other. As a re- 
sult, these bands are rather flat, disentangled and only 
crossing on limited symmetry points or lines. We can find 
a pseudo-gap <~ 0.2 eV in the vicinity the Fermi level, 
which is only closed around the L (5,5,5) point on the 
trigonal axis of the COO-I state. 

Figure 5 shows the band structure and the density 
of states for the COO-II state. The band structure of 
the COO-II state is similar to that of the COO-I state 
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Fig. 4. Band structure and the density of states for the COO-I 
state obtained with V = 0.1 eV, t pd = 0.35 eV, D = 0.25 eV and 
ctg = 1, !)■ The broken lines in the figure denote the Fermi 

energy. The symmetry points K(|,0, f), L(|,|,|), U(±,±,l), 
W(f ,0, 1) and X(0,0,1) arc chosen. 



(see Fig. 4), and band gaps at the L points which are 
on the direction from the centre of the Fe 4 04 cube to 
the electron-rich sites are small; for instance, gaps at 
the L(i, 5, 5) and L(— |, —5, 5) points are small for the 
COO-II state depicted in Fig. 3(b). Again, a pseudo-gap 
in the vicinity of the Fermi level can be found among 
the low-lying four bands. The gap is closed only on the 
r— K symmetry line along the [110] c direction, which is 
on the mirror plane of the COO-II state, and thus the 
COO-II state is a zero gap semiconductor. The density 
of states near the Fermi level reduces with increasing V. 
Since these COO states are on the verge of gap opening, 
it is expected that further symmetry lowering caused by 
lattice distortion in the low-temperature phase would re- 
sult in the formation of the insulating gap. In § 5, we will 
discuss the electronic state in the low-temperature phase 
in relation to the lattice deformations. 
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Fig. 5. The same as Fig. 4 but for the COO-II state obtained with 
«S = 75(M,0). The symmetry points K(f,f,0), L(±,i,i), 
U(|,l, i) W(i,l,0) and X(0,1,0) arc shown. 



4. Formation Mechanism of COO State 

4-1 Interacting Fa tetrahedral molecules of B sites 

In the previous section, we have found the COO states 
stabilised in large t p d regime in the spinless three-band 
Hubbard model on the pyrochlore lattice. However, it 
is not obvious why such unconventional orbital order- 
ings realise in this system. To gain more insight about 
these peculiar orbital orderings, let us consider an iso- 
lated tetrahedron of B sites in infinite U limit as the 
first step. The pyrochlore lattice can then be treated as 
a fee lattice of these tetrahedral molecules of B sites. 

Within HF approximation, only one kind of orbital can 
be occupied in each of the B sites, if U = 00. We fur- 
ther assume the solution possesses the trigonal symmetry 
with the [lll]c axis and is invariant with combined oper- 
ator of time reversal and vertical mirror. The symmetry 
is the same as the COO-I state and here, we use the same 
site labels in Fig. 3(a). With these assumptions, the wave 
function at site 1 can take one of the a\ or e ± symmetry 
orbitals: 



|l,e+> = ^= 



II, e" 



I 
I 

V3 
-1 

7! 



(\l,yz) + \l,zx) + \l,xy)) 
(uj 2 \l,yz) +w\l,zx) + \l,xy)) 
{u\l,yz} +w 2 \l,zx) + \l,xy)) 



and the wave function at site 2, 3 and 4 can be written 

as 

1 2) = v\2,yz) + u\2,zx) + u*\2,xy) 

1 3) = u* |3, yz) + v\3, zx) + u\3, xy) 

|4) = u\A,yz) + u*\A, zx) +v\A,xy), 

where u and v are complex and real constants, respec- 
tively, and satisfy the relation 2\u\ 2 + v 2 = 1. The ai and 
e ± symmetry molecular orbitals (MO) on the triangular 



cluster of sites 2, 3 and 4 are described as 
|234,oi) = -^(|2) + |3) + |4)) 

|234, e +) = -J=(u, 2 |2)+u;|3> + |4)) 

|234,e-> = ^(^|2) +^ 2 |3> + |4)) 

and if we express one-body energies for these a\ , e + and 
e~ symmetry MOs as e°, e +1 and e _1 , respectively, e™ 
can be described as 

s n = 2Re[w™(3|iJ|2)] 

= 2Re[co n {-(u*)% + u 2 t pd + v 2 t pd }}, 

where H is the Hamiltonian described in eq. (1) and t a = 
3|idd<r|/4. Here, for simplicity, the values of parameters 
tdd-m D, C and V are assumed to be zero. Substituting 
u = e 1 ^ cos#/\/2 and v — sin 9, we obtain 

e 11 = — t a cos 2 9 cos \2(j)~ ~^-n 



2tt 



+ tpd cos 9 cos I 2<j) 



+ 2t p d sin 2 9 cos ( 



(4) 



Either e + or e symmetry MO can have the lowest en- 
ergy among the three MOs on the triangular cluster: 



£ ±X = -yjtl + tvtpd + t 



2 

pd 



at 9 — and 



2(f> = ± 



tan 



2tpd t(j 3 



The value of 2(f) for this lowest energy MO is within it/ 3 < 
1 20| < 27r/3 and thus orbitals on each sites are complex 
orbitals. For example, the orbitals at sites 2, 3 and 4 are 
x + , y + and z + orbitals, respectively, when 9 = and 
2(f> — ir/2 are chosen for the e + symmetry MO and this 
orbital alignment on the three sites is the same as that for 
the COO-I state in Table I. These MOs on the triangular 
cluster are further hybridised with the orbitals on site 1 
having the same symmetry. The hopping integrals, i.e., 
h° = (l,ai|^|234,ai) and h* 1 = (1, e ± |i?|234, e ± ), can 
be written as 

(2-7T 
— 3~ ? 

and thus the energies of the bonding e£ and antibond- 
ing states of the MO and orbital on site 1 in each 
symmetries can be obtained as 

1 



£ b,a — 



e n ± V(e") 2 +4(/i n ) 2 



(5) 



For the tetrahedron with two electrons, the lowest- 
energy electronic configurations are those with an elec- 
tron on the e^ symmetry MO of the triangular portion of 
the tetrahedron (E ) and the other electron on the bond- 
ing MO with e T symmetry (E^) , i.e., E+E~ or E~E£ 
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Table II. The largest four eigenvalues Wi and the symmetries 7; 
and parameter values <pi and \i) of corresponding eigen- 
functions of the density matrix of a tctrahcdral molecule in the 
COO-I state on the pyrochlore lattice within the HF approxima- 
tion. The energies of Bloch states Si with 7, symmetries at the T 
point are also shown. In (a), {7=40.0 eV and t^d-n = D = eV, 
and in (b), (7=4.0 eV and D = 0.25 eV are adopted; t dp =0.4 eV, 
V = eV and the other parameters are the same as § 3. 





It 


Wi 


0i (°) 


2& (°) 


Xr (°) 


Si (eV) 


(a) 


E+ 


0.806 


-11.7 


93.3 


0.3 


-1.032 






0.768 


-12.9 


95.0 


-58.6 


-0.721 




A 1 


0.382 


-11.5 


95.0 


0.7 


0.008 




E- 


0.044 


-13.1 


92.6 


31.3 


1.760 


(b) 


E+ 


0.800 


-10.2 


91.3 


4.2 


-0.791 






0.727 


-23.2 


112.8 


-54.9 


-0.637 




A l 


0.414 


-11.2 


89.4 


4.1 


-0.007 






0.036 


-19.9 


91.9 


34.7 


1.061 



configurations. Note that configurations such as E^E^ 
cost large Coulomb repulsion energy at site 1 and are not 
favourable. The minimum total energy of E + E^ config- 
uration with respect to 9 and (j> is found at 9 <~ —60° and 
2cf) ~ 89° with parameter values tdda = —0.275 eV and 
t p d = 0.4 eV, where the COO-I state is the ground state 
on the pyrochlore lattice. 

Now let us consider electron hopping between tetra- 
hedral molecules on the pyrochlore lattice in addition to 
the internal hopping in each of them. The Hamiltonian 
for the Bloch function with wave vector k can be easily 
obtained by replacing the hopping matrix of the isolated 
tetrahedron: 



where site f is equivalent site to j on the adjacent tetra- 
hedron sharing a (110) c chain with site j and (ry) 
represents the position of the site i (j'). For the Bloch 
states with k around the T point, the intra- and inter- 
molecular hopping terms are in-phase (~ 2£^ I/ ), and thus 
the energy gain from the hopping term is approximately 
two time as large as that of the isolated tetrahedron. 
This stabilises a COO state, where the lowest two MOs 
in each tetrahedron arc mainly occupied, and they are co- 
herently aligned. On the other hand, for the Bloch states 
with k in the vicinity of the boundary of the first Bril- 
louin zone, the phases of the intra- and inter-molecular 
hopping terms are not matched. This causes mixture of 
the lowest two MOs and the higher ones through the 
inter-molecular hopping, because of breaking of the local 
symmetry. This would modify the ordered MOs of the 
COO state from those in the isolated tetrahedron. 

To verify the above scenario with interacting tetra- 
hcdral molecules of B sites, the density matrix Pi^j v = 
(Cj U c in) f° r a tetrahedral molecule of B sites (in a 12 x 12 
size) in the COO-I state on the pyrochlore lattice is cal- 
culated and diagonaliscd. Results obtained within the 
HF approximation with two kinds of parameter sets are 
shown in Table II: in (a), [7=40.0 eV and t dd7T = D = 0, 
and in (b), [7=4.0 eV and D = 0.25 eV; in both calcu- 



lations, tpd — 0.4 eV and V = eV are adopted and the 
other parameters are the same as the previous section. 
The largest four eigenvalues Wi and the symmetries 7, 
and parameters (0,, 4n and Xi) °f corresponding eigen- 
functions are listed. Here, Xi is the mixing angle between 
the wave functions at site 1, and the MO on the trian- 
gular part is defined as 

|7i) = sin x» 1 1,7*) +cosxi|234,7i). 

The energies of Bloch states Ei with symmetry 7^ at the 
T point are also shown. 

The values of 9 and 4> are almost the same among 
the eigenfunctions in Table 11(a). This is because only 
one particular orbital can be occupied in each sites in 
large U limit. Variations in these parameters are larger 
for the results with realistic parameters in Table 11(b). 
However, no qualitative differences are found between 
the two. As expected, the lowest-energy two MOs E + 
and E^ are dominantly occupied. In addition, about 0.4 
electron resides on the Ai symmetry MO with the third 
lowest energy. This is because the lowest two MOs can 
be hybridised with the higher ones through the inter- 
molecular hopping. The value of 9 ~ —12° is much 
smaller than that in isolated tetrahedron: 9 ~ —60°. 
This large discrepancy in the value of 9 can be ex- 
plained by the differences in the MO occupations. In- 
deed, if the lowest-energy two MOs and the third-lowest 
MO are assumed to accommodate 0.8 and 0.4 electron, 
respectively, in the isolated tetrahedron, the configura- 
tion (E + ) a - s (E h )°- s (A 1 ) 0A has the minimum total en- 
ergy at 9 <~ —11° and 2(f> ~ 94° with parameter values 
tdda = -0.275 cV and t pd = 0.4 eV. 

We have described the COO-I state as the ferro-MO 
ordered state of the tetrahedral molecules of the B sites, 
assuming the trigonal symmetry. Since there are similar- 
ity among the COO states, one may consider all these 
states can be described by a more general scheme. In 
the next subsection, we will find that all these COO 
states can be regarded as a ferromagnetic state of fic- 
titious orbital moments of the pseudo-<i symmetry MOs 
on the tetrahedral molecules. In this context, the differ- 
ences among them are only magnetisation directions of 
the ferromagnetic state. 

4-2 Stability of COO state and pseudo-SO(3) symmetry 
A COO state in the e g orbitals has been proposed to 
explain the properties of ferromagnetic metallic phase 
of doped manganites, where the static Jahn- Teller de- 
formation disappears. 30-32 The e g orbitals \x 2 — y 2 ) and 
|3z 2 — r 2 ) are combined with the complex number co- 
efficients, e.g., ^ [|3z 2 - r 2 ) ± i|a; 2 - y 2 )] in this COO 
state, and because of the isotropic charge distribution, 
the COO state is not associated with the Jahn- Teller de- 
formation. Kubo and Hirashima, however, have argued 
that the COO state in ferromagnetic metallic mangan- 
ites is unstable, when the electron correlation effect is 
taken into account. 33 We, therefore, performed calcula- 
tions on a finite size system using the exact diagonali- 
sation method to clarify whether the COO states in our 
three-band model are stabile under the electron corrcla- 
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Fig. 6. Orbital orders found in the ground state of the 16- 
site model obtained with the exact diagonalisation method are 
mapped on the V-t p d plane. In the calculations, D = 0.25 eV 
and cng = "t|(1: b 1) are assumed. 



tion effect and why these pecufiar COO states are pre- 
ferred. 

The three-band model described in § 2 is employed 
for the finite size system, where sixteen B sites in the 
conventional fee unit cell accommodating eight electrons 
are considered with the periodic condition. To reduce the 
number of the many-body basis functions, we assume 
U = oo. The other parameters are the same as § 3. The 
Hamiltonian of the 16-site model is diagonalised with 
the Lanczos algorithm 34 to obtain the ground-state wave 
function. 

The orbital orders found in the ground state of the 
16-site model are mapped on the V-t p d plane in Fig. 6. 
In the calculation, D = 0.25 eV and as = ^= (1, 1, 1) are 
adopted. The regions of the ground states corresponding 
to the ROO-I and COO-I states can be seen. To examine 
instability toward the Verwey charge ordering, the stag- 
gered potential Y] ifl Si was also applied in the ROO-I 
region, where Si = +S and Si = —6 for the [110] c and 
[110] c chains, respectively, and 2<5=0.05 eV is adopted. 
The difference between the occupation numbers on the 
[110] c and [110] c chains continuously varies from ~ 0.1 to 
~ 0.3 with increasing values of V and t p d- These results 
are consistent with the phase diagrams in Fig. 2 obtained 
from the HF calculations. 

As was discussed in § 3, the energies of the three COO 
states are nearly degenerated within ~0.01 eV, and each 
of them can be stabilised when as is oriented to the di- 
rection of their total magnetic moments. In fact, almost 
the same diagrams as in Fig. 6 were found except the 
COO-I region is replaced by those for the COO-II or 
COO-III states with a s = 1,0) or a s = (0,0,-1), 

respectively. These COO states can even continuously 
transform into the others as we gradually vary as- As 
discussed later, the COO states are results of sponta- 
neous breaking of the inherent pseudo-spherical symme- 
try analogous to the ferromagnetic order in the Heisen- 
berg model. 

Figure 7 shows the relationship between the trigonal 



t > 1 > 1 > r 




I i i i i i i i I 

0.1 0.2 0.3 

D(eV) 

Fig. 7. Trigonal field D dependence of the total orbital moment 
(fiB per B site) obtained with V = 0.1 eV, t p d = 0.35 eV and 

«s = 5i(i,i,i). 



field D and the total orbital moment per B site. While 
large orbital moments are induced and the COO-I state 
is found to be the ground state for D > 0.1 eV, the COO 
state collapses and only small orbital moments appear for 
small values of D. From the analysis of the 3x3 density 
matrix {c\ v c itl ) on each site, we found orbital occupation 
in each site for D < 0.1 eV is described as the mixed 
state, where two kinds of the e£ orbitals are occupied 
with the same probability, except for a small deviation 
(<~ 10%) caused by the spin-orbit interaction. Thus, this 
ground state can be considered as an orbital-liquid state 
with small orbital polarisation owing to the spin-orbit 
interaction. On the other hand, for D > 0.1 eV, one par- 
ticular orbital shown in Table I is dominantly occupied in 
each site, and the occupation probabilities for the other 
two are small (less than 5% for D = 0.25 eV). Contrast- 
ingly, the COO-I state remains to be the lowest energy 
solution with the HF approximation down to D=0. It 
is also noteworthy that even when the spin-orbit cou- 
pling £ on all site is set to zero except for a site and as 
is oriented to the trigonal axis of this site, the COO-I 
state is still obtained with D > 0.1 eV. This also shows 
robustness of the COO-I state for D > 0.1 eV. 

To understand this change in the ground state as a 
function of D, let us consider an isolated tetrahedral 
molecule of the B sites accommodating one electron. The 
tig orbitals on four B sites in the tetrahedron are reduced 
to five different basis sets of irreducible representations 
in the T d symmetry: A x + E + T x + 2T 2 . The MOs arc 
classified according to the symbols of the irreducible rep- 
resentations, and the bases and energies are described in 
the Appendix. The energy levels corresponding to the 
bonding and antibonding states of two T 2 basis sets (T!f 
and T 2 ) are labelled T 2 t and T 2a , respectively. In Fig. 8, 
MO level diagrams for the isolated tetrahedral molecule 
obtained with t p d — eV (a) and t p d — 0.3 eV (b) are 
shown. 

For small value of t p d, the lowest is the A\ level, where 
the a\ g orbitals are occupied by 0.25 electron in each. 
The second lowest is the E level for small values of D 
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Fig. 8. Molecular-orbital levels of the isolated tetrahcdral 
molecule of the B sites as a function of D obtained with t pl i=0 eV 
(a) and i p( j=0.3 eV (b). 



and is the T 2b level for D > 0.45 eV. For large value of 
t p d, the occupation of the orbitals is more favourable, 
and the E level, where the e£ orbitals are occupied for 
all B sites, is the lowest level. However, the T 2 fc level 
becomes the lowest for D > 0.45 eV. 

From these level diagrams, we can explain the change 
in the orbital ordering qualitatively. For the isolated 
tetrahedron with two electrons, the E 2 closed-shell con- 
figuration is favourable in the case that the energy sep- 
aration between the lowest E and higher levels is large. 
This costs large on-site Coulomb energy since two kind 
of the orbitals, i.e., e u and e„, should be occupied 
on each B sites (see Table A-l). In this situation, large 
quantum fluctuation is expected even on the pyrochlorc 
lattice, since electrons should avoid occupying the same 
site with this limited orbital degree of the freedom. On 
the other hand, if the T 2 b and E levels are nearly degen- 
erated {e(T 2 b)—e(E) < 0.2 eV], electrons would prefer to 
occupy a particular orbital in each of the B sites utilis- 
ing both the E and T 2 b basis sets to minimise the on-site 
Coulomb energy. 

From above discussions, the pseudo-degeneracy of the 
E and T 2 b levels is the requisite for the formation of 
the COO state. If we regard these MOs as irreducible 
representations of a d state in the T<j symmetry, i.e., d = 
E + T 2 b, the fictitious d symmetry MOs with magnetic 
quantum number m can be defined as 

\d,m = ±2) = ^=(\E,v)±i\T 2b ,Q), 



\d,m- 



V2 

-±1) = ^=(T\T 2b , V ) 



i|22t,0). 



\d,m = 0) = \E, u). 



(6) 



If the pseudo-S'0(3j symmetry in the Coulomb interac- 
tions among these d MOs are further assumed, the most 
lowest energy configuration for accommodating two elec- 
trons in this spinless d MOs is that obeys Hund's second 
rule: the d 2 configuration, where \d,m = 2) and \d,m = l) 
orbitals occupied with the fictitious total orbital mo- 
ments L = 3. Thus, the COO state can be interpreted as 
a ferromagnetic (FM) alignment of the fictitious orbital 



moments of the d MOs, which is parallel to the real mag- 
netisation direction of the COO states. From this view- 
point, the differences among the COO states are merely 
the directions of the magnetisation of the FM state. The 
COO-I, COO-II and COO-III are, therefore, only three 
particular cases of the FM state with magnetisation di- 
rections (lll)c, (110)c and (001)c, respectively. 

To check the validity of the assumption of the pseudo- 
5*0(3) symmetry, the density matrix Pi^j v for the tetra- 
hedral molecule is diagonalised. For the COO-III state 
corresponding to the FM state with the [001] c mag- 
netisation direction, the largest three eigenvalues are 
wi = 0.73, w 2 = 0.54 and w 3 = 0.40. The overlap 
between these eigenfunctions of p and the d orbitals in 
cq. (6) are | (wi \d, m = 2) | = 0.96, | (w 2 \d, m = 1) | = 0.88 
and | (103 1 d, m = 0)| = 0.51. As expected, the dominantly 
occupied two MOs are correspond to the \d, m = 2) and 
\d,m = 1) orbitals and their deviations from these d or- 
bitals are fairly small. The COO-II state correspond to 
the FM state with the [110] c direction. The wave func- 
tions of the d states can be easily obtained by rotating 
the quantisation axis to the magnetisation direction us- 
ing the representation matrix D m , m i [R(a, (3, 7)] of the 
rotation R(a, (3, 7) described by the Eulerian angles. For 
instance, the \d,m') orbital with the [110] c quantisation 
axis can be obtained as 
2 

\d,m' = 2) = \d,m) D m!2 [R(Tr/4,TT/2,0)] 



-^(\T2b,0 



\T2b,v)) 



1 



\d,m' = l) = -j=\E,v)+ 1 - 



2^2 

T 2b ,0 + \T 2b ,ri)) 



\T2b,Q, 



\d,m' = 0) = ~\E,u) + ^\T 2b ,0. 

The largest three eigenvalues of p in the COO-II state are 
w\ = 0.78, w 2 — 0.66 and W3 = 0.39 and the overlaps 
are \(wi\d,m' = 2)| = 0.99, \(w 2 \d,m' = 1)| = 0.87 and 
I (u>3 1 d,m' = 0)| = 0.39. Again, dominantly occupied two 
MOs are reproduced well with the pseudo-5 , 0(3) sym- 
metry assumption. The examination of p in the COO-I 
state leads to the same conclusion, and we now realise 
that the E + and E^ MOs discussed in the previous sub- 
section correspond to the \d,m = 2) and \d, m = l) MOs, 
respectively. 

The pseudo-spherical symmetry, which permits the 
creation of Hund's rule state on each Fe4 unit, is the fun- 
damental ingredient for the formation and stabilisation of 
the COO states in the present system. In contrast, such 
inherent symmetry has not been discussed and seems to 
be absent in the e g system for the ferromagnetic man- 
ganites. 30-33 

5. COO State in the Low- Temperature Phase 

5. 1 Lattice distortion and related properties in the low- 
temperature phase 
In the low-temperature phase (T < Ty), the lattice 
distortion to the monoclinic (Cc) or triclinic (PI) symme- 
tries with a cell size V2a c x V2a c x 2a c occurs. However, 
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Fig. 9. Illustration of an a c /\/2 X a c /\/2 X 2a c cell for the Pmca 
and Pmc 2i pseudo-orthorhombic structures. 4, 35 The circles rep- 
resent Fe ions on the B sites. The six inequivalent B sites in the 
Pmca structural model arc labelled Bla, Sib, B2a, B2b, B3 and 
BA after ref. 4, and those in the parentheses are for the Pmc 2\ 
structural model in ref. 35. Prominently contracted Fc-Fc bonds 
~ 0.1 A in the Pmca (those with the black and grey arrows) and 
in the Pmc2i (those with the black arrows) structural models 
arc also shown. The monoclinic a m -, fe m - and along 
the [110]c, [110] c and [001] c directions, respectively. 



attempts to determine the true crystal structure have 
been hindered by overlapping of the Bragg peaks owing 
to the presence of pseudo-symmetries. As an approxi- 
mated crystal-structure model for the low-temperature 
phase, Iizumi et al. proposed an a c /v / 2xa c /v / 2x2a c sub- 
cell imposing orthorhombic symmetry constraints on the 
atomic positions to reduce the number of free parameters 
from the true \f2a c x \pla c x 2a c Cc supercell accommo- 
dating 16 inequivalent B sites in their neutron diffrac- 
tion study. 35 They used the centric Pmca and acen- 
tric Pmc2i space groups for candidates of orthorhom- 
bic pseudo-symmetry. Recently, Wright et al. refined the 
crystal structure of the low-temperature phase from x- 
ray and neutron diffraction experiments. 3,4 The same 
crystal structure model with the orthorhombic Pmca 
pseudo-symmetry is applied for the analysis. The in- 
equivalent sites are reduced from 16 sites to 6 sites within 
the approximation, and they are labelled £?la, Bib, B2&, 
B2b, B3 and B4 (see Fig. 9). Jahn- Teller-like distortions 



of the oxygen octahedra surrounding the Bl sites are 
found below Ty from their experiments, and charge dis- 
proportionation is inferred from variations in the Fe-0 
bond lengths. 

The electronic state in the low-temperature phase has 
been investigated using LDA + U calculations 16-19 based 
on the monoclinic lattice structure shown by Wright et 
al. 4 From the LDA + U calculations, a state with simul- 
taneous charge and antiferro-orbital ordering is obtained. 
In the state, the £?2a, B2b and J53 sites are vacant, and 
the xy orbitals are occupied at the B4 sites. The yz and 
zx orbitals align alternately along the Bl chains, where 
the Jahn- Teller-like distortions are observed. 

The cell distortion of magnetite in the low-temperature 
phase is closely approximated by that with a rhombo- 
hcdral symmetry, 2,36,37 although the atomic displace- 
ments in the cell have pseudo-orthorhombic symmetry. 4 
The monoclinic distortion in the low-temperature phase 
has been conventionally interpreted as the result of the 
combination of the microscopic orthorhombic distortion 
caused by charge ordering and macroscopic rhombohc- 
dral distortion owing to the magnetostriction. The mag- 
netic easy axes are oriented along the [11 l] c direction 
in the high-temperature phase and the monoclinic c m - 
axis (more precisely, it tilts toward the [Hl] c direction 
by 2° 38 ) in the low-temperature phase. 39 Here the mon- 
oclinic a m -, b m - and c m -axes are along the [110] c , [110] c 
and [001] c directions, respectively. Interestingly, the easy 
axis switching accompanied by simultaneous change of 
c m -axis is observed by applying external magnetic field 
below Ty, 40 and the dependence of the Verwey transition 
temperature under the application of magnetic field was 
also reported. 41 ' 42 Moreover, the magnetoelcctric effects 
which indicate triclinic distortion in the low-temperature 
phase were also discussed. 43-48 

As was discussed in the previous section, the COO 
state is a kind of magnetic order within the orbital de- 
gree of freedom. This means that its orbital polarisa- 
tion depends on the magnetisation directions and thus 
strong coupling between the magnetism and lattice dis- 
tortions, i.e., magnetostriction, is anticipated. Because 
of this strong coupling with the lattice, the COO state 
in the high-temperature phase would be modified by the 
monoclinic lattice distortion below Ty. In the remainder 
of this section, we will discuss the relationship between 
the COO state and lattice distortion. A COO state ex- 
pected to be stabilised by the monoclinic lattice distor- 
tion in the low-temperature phase will be also presented. 

5.2 Relationship between lattice distortions and COO 
states 

To discuss the relationship between the lattice defor- 
mation and COO states, the hopping integrals between 
the dominantly occupied orbitals hij = (i, n\H t \j, v) in 
the COO-I, COO-II and COO-III states are examined. 
Here, \i,fi) represents the dominantly occupied orbital \x 
at site i which is listed in Tabic I. If the absolute value 
of the hopping integral is increased by the Fe-Fe bond 
contraction, the energy of corresponding bonding state 
is decreased, and this would reinforce the stability of the 
ordering. In Table III, \hij\ expressed in terms of t a and 
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Table III. Expressions in terms of t a and t p d and the first 
derivative with respect to t„ of the absolute value of the hop- 
ping integrals between the dominantly occupied orbitals h{j = 
{i,H\H t \j,v) in the COO-I, COO-II and COO-III states. 
denotes the dominantly occupied orbital /i at site i in the nota- 
tion in Fig. 3. Here, t a = 31^0-1/4, and t^-n la neglected. Ht is 
the hopping term in eq. (1). 



state 


hij 


1 hij 1 


d\hij\ 
dt„ 


COO-I 


(3,y+\H t \2,x+) 
(4,z+\H t \l,e-) 


\J 6 tpd 


0.50 



COO-II 


(3,110|i/ t |2,110) 
(3,110|5 t |4,z+) 


\{2t a + tpd) 

-^(ta + 2t pd ) 


0.67 
0.29 




(4,z+\H t \l,z-) 


tpd 





COO-III 


(2,lll|flt|l,lll) 
(4,lll|5 t |l,lll) 


\{t a + 2t pd ) 
3 (t(T tpd) 


0.33 
0.33 



t p d, and the first derivative of \hij\ with respect to t a are 
listed. Here, t a — 31^0-1/4 is the absolute value of the 
hopping integral of the a-bond, which is most sensitive to 
the change in the bond length among the hopping inte- 
grals. Note that there are the relations /i32=/i24=^43 and 
hi4=h\2=hi3 for the COO-I state, /i34=/i42=/i2i=^i3 for 
the COO-II and COO-III states, and hi 4 =h 2 3 for the 
COO-III state. 

Judging from the derivatives, the Fe-Fe bond between 
the two electron-poor sites of the COO-II state, on which 
orbitals 1 110) are occupied, is most sensitive to the vari- 
ation of t a . Thus, the orthorhombic distortion, in which 
the Fe-Fe bond along the [110] direction is contracted, 
is expected to stabilise the COO-II state. Similarly, the 
COO-I state is accompanied by the rhombohedral lat- 
tice distortion; three Fe-Fe bonds on the triangle of sites 
2, 3 and 4 of each of the Fe404 cubes are compressed, 
and this leads to elongation of the cubic cell along the 
[lll] c diagonal direction. On the other hand, only small 
tetrahedral distortion is expected for the COO-III state. 

5.3 Hartree-Fock calculations for the low-temperature 
phase 

The spinless three-band Hubbard model described in 
cq. (1) was investigated to find the electronic state re- 
alised in magnetite below Ty. The calculations were car- 
ried out within the HF approximation [eq. (3)]. In order 
to introduce the effects of the lattice deformation for the 
low-temperature phase, we considered the variation in 
the values of hopping integrals as a function of Fe-Fe 
bond length r. For simplicity, we take into account such 
r dependence on the value of the hopping integral for 
the cr-bond — 1„ = Stddo/^, which is most sensitive to 
the change in r among the hopping integrals. The value 
of hopping integral is approximately proportional to r~ 5 , 
and we expand t a up to the first order of the deviation Sr 
from the Fe-Fe distance in the high-temperature phase 
r c = 2.966 A as 



Table IV. Deviation of the bond length from that in the high- 
temperature phase Sr and corresponding value of the hopping 
integral t a (= 31^^3-1/4) adopted in the HF calculations are 
shown for each of the Fe-Fe bonds among the B sites. Those 
evaluated from two sets of atomic coordinates in the Pmca 4 ' and 
Pmc 2i 35 structural models with /3 CT = 3 are presented. 



bond 


Pmca 
Sr (A) 


U (eV) 


bond 


Pmc2\ 
Sr (A) 


to (eV) 


51- 


-51 


0.000 


0.206 


51- 


-51 


-0.046 


0.255 


51- 


-53 


-0.107 


0.318 


51- 


-51 


0.046 


0.158 


51- 


-54 


-0.013 


0.220 


51- 


-52 


-0.050 


0.259 


52- 


-52 


0.000 


0.206 


51- 


-53 


0.052 


0.152 


52- 


-53 


0.080 


0.123 


51- 


-55 


-0.118 


0.330 


52 


-54 


0.042 


0.163 


51 


-56 


-0.008 


0.215 


53 


-54 


-0.082 


0.292 


52- 


-53 


-0.002 


0.209 


53- 


-54 


0.083 


0.120 


52 


-53 


0.005 


0.201 










52- 


-54 


-0.004 


0.211 










53- 


-54 


-0.065 


0.274 










54- 


-54 


-0.079 


0.289 










54- 


-54 


0.079 


0.123 










54- 


-55 


0.117 


0.085 










54- 


-56 


0.080 


0.123 










55- 


-56 


-0.104 


0.315 










55- 


-56 


0.105 


0.097 



t a = £(1 - 5p a Sr/r c 



where t% = 0.206 eV is the hopping integral for the high- 
temperature phase. A constant parameter f3 a is intro- 
duced to change the extent of the distortion effects for 
the cr-bonds, and we assumed f3 a = 3 in the calculations. 

The a c j\[2 x a c j\[2 x 2a c orthorhombic unit cell con- 
sists of 16 5 sites (see Fig. 9) is assumed. For simplicity, 
effects of atomic displacements in the cell are considered 
only as variations in hopping integrals, and the mono- 
clinic tilt of the c m -axis toward the a m -axis are neglected. 
Note that the pairs of sites labelled with 'a' and 'b', e.g., 
51a and 51b, arc equivalent without this monoclinic tilt 
of the c m -axis because of recovery of the 6 m -c m mirror 
plane. Two sets of the atomic coordinates with the cen- 
tric Pmca 4 and acentric Pmc 2\ 35 symmetries were used 
to evaluate the values of t a , and those adopted in the HF 
calculations are listed in Table IV. The direction of the 
minority spin moment is assumed to be along the c m - 
axis, i.e., £*s = (0, 0, —1) since the magnetic easy axis for 
T < Ty is approximately along this direction. The other 
parameters are the same as those used in the calculations 
for the high-temperature phase: D = 0.25 eV, U — 4.0 eV 
and C = 0.052 eV. The wave vectors throughout the first 
Brillouin zone of the orthorhombic lattice were sampled 
in the HF calculations. 

From the both calculations using the hopping integrals 
evaluated from the Pmca and Pmc2\ structural mod- 
els, we obtained a charge and complex-orbital ordered 
state (COO-IV). Schematic representation of the COO- 
IV state is shown in Fig. 10, and its details are listed in 
Table V. For the both symmetries, the COO-IV state is 
found as the lowest energy solution in almost the same 
area on V-t p d plane where the COO phases presence in 
Fig. 2(c). Although we assumed the orthorhombic lattice 
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Fig. 10. (a) Schematic representation of the COO-IV state found in the HF calculations for the low-temperature phase. The large and 
small circles denote the Fe ions on the B sites and oxygen ions, respectively. The black and grey arrows signify the directions of orbital 
moments with large and small values of occupation number at the B sites, respectively, (b) The COO-IV state projected on the [lll] c 
plane. The circles represent the Fe ions on the B sites. The open and hatched circles inside the triangles are above and below kagome 
plane, respectively, and the others are on the kagome lattice. The site labels in the Pmca symmetry are also shown. 



structure, the symmetry of the COO-IV state is sponta- 
neously lowered from the Pmca and Pmc2\ to the mono- 
clinic P2/c and Pc, respectively. However, the symmetry 
of the COO-IV state is essentially P2/c and deviation 
from this is small in the calculation with the Pmc2\ lat- 
tice symmetry. Thereby, for the remainder of this paper, 
we refer the B sites by the labels for the Pmca lattice 
symmetry 

As seen in Fig. 10, the COO-IV state can be regarded 
as an alternative stacking of microscopic layers (thickness 
a c ) of the COO-II states with two different orientations: 
those having easy-axes along the [101] c and [011] c direc- 
tions. The horizontal planes at the z-coordinate of the 
lattice z m = 1/4 and z m = 3/4 correspond to 'domain 
walls', on which the B2 sites are located. Each of the B2& 
sites is placed at the intersection of a [011] c chain of the 
1 1 1 ) orbitals and a [101] c chain of the 1 101) orbitals, and 
its orbital moment is along the [1 1 l] c direction. Similarly, 
each of the B2b sites is positioned at a intersection of a 
[011] c chain of the Ix*) orbital and a [101] c chain of the 
jy*) orbitals, and its orbital moment is oriented along 



the [110] c direction. The B1-B3 bonds are the most con- 
tracted bonds, and both the chains of the 1 1 1) and 1 101) 
orbitals of the COO-II layers are on these bonds. From 
the discussions in § 5.2, this explains why the COO-IV 
state is stabilised with the lattice distortion in the low- 
temperature phase. 

The magnetic easy axis of the COO-IV is oriented to 
the [223] c direction tilted from the c m -axis toward the 
[lll] c direction about 43° and the total orbital moment 
<~ 0.2 per B site. The experimental easy axis is also 
tilted from the c c -axis toward the same direction below 
Ty. However, the tilting angle is about 2°. 38 

The charge disproportionation in the COO-IV state 
is rather moderate similar to the COO-II state, and the 
charge order satisfies the Anderson condition, which con- 
strains every tetrahedron to have two electrons. From 
the bond valence sums analysis of the x-ray and neutron 
diffraction experiments, 3 ' charge ordering with appar- 
ent local charges of Fe 2 4+ and Fe 2 6+ has been found in 
the low-temperature phase. Although the small magni- 
tude of charge splitting is consistent with our results, the 
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Table V. Details of the COO-IV state obtained from the calculations with the hopping integrals evaluated from the Pmca and Pmc2\ 
structural models. The z-coordinate of the atomic position (z m ), dominantly occupied orbital and its orbital moment (in units of /is), 
and occupation numbers (for V = 0.3 eV and t p[ j = 0.35 eV) on each of the 16 sites are listed. 





orbital 


LJ 


= cxp(27ri/3) ] 


moment 


(n) (Pmca) 


(n) (F 


mc2 1 ) 





1 101) = 


V6 


2\zx) -i\xy) + i\yz)] 


|(1,0,1) 


0.411 


(Bib) 


0.418 


(Bib) 





\y-) = 


1 

V2 


\xy)+i\yz)} 


(0,-1,0) 


0.617 


(Bid) 


0.612 


(Bla) 


i 

8 


1 101) = 


7E 


2\zx) -i\xy) +i\yz)] 


1(1,0,1) 


0.434 


(S3) 


0.463 


(B2) 


1 

8 


\y + ) = 


i 


\xy) - i\yz}] 


(0,1,0) 


0.539 


(B4) 


0.519 


(S3) 


1 

4 


|iii) = 


1 


\xy) + u) 2 \yz) + lv\zx)] 




0.486 


(B2a) 


0.453 


(B4a) 


1 

4 


1 110) = 


1 

V& 


2\xy) + i\yz) - i\zx)] 


-|(1,1,0) 


0.540 


(B2b) 


0.570 


(B4b) 


3 

8 


|011) = 


V6 


2\yz) + i\xy) -i\zx)] 


1(0,1,1) 


0.434 


(S3) 


0.432 


(£5) 


3 
8 


\x+) = 


1 

V2 


\zx) - i\xy)] 


(1,0,0) 


0.539 


(B4) 


0.533 


(B6) 


1 

2 


\x~) = 


1 
V2 


]zx) + i\xy}] 


(-1,0,0) 


0.617 


(Bla.) 


0.612 


(Bla,) 


2 


|011) = 


75 


2\yz) + i\xy) - i\zx}] 


1(0,1,1) 


0.411 


(Bib) 


0.418 


(Bib) 


5 
8 


|011) = 


i 

V6 


2\yz) + i\xy) - i\zx}] 


§(0,1,1) 


0.434 


(B3) 


0.463 


(B2) 


5 
8 


\x+) = 


1 

V2 


\zx) - i\xy)] 


(1,0,0) 


0.539 


(B4) 


0.519 


(S3) 


3 
4 


|111) = 


1 

V3 


\xy) + Lo 2 \yz) + u\zx)} 


A* 1 ' 1 ' 1 ) 


0.486 


(B2a) 


0.453 


(54a) 


3 
4 


1 110) = 


1 

>/5 


2\xy) +i\yz) -\\zx)\ 


-§(1,1,0) 


0.540 


(B2b) 


0.570 


(B4b) 


7 

8 


\y + ) = 


1 

V2 


\xy) - i\yz)] 


(0,1,0) 


0.539 


(B4) 


0.533 


(B6) 


7 

8 


1 101) = 


1 

V5 


2\zx) - i\xy) + i\yz}] 


1(1,0,1) 


0.434 


(S3) 


0.432 


(B5) 



charge order found in the analysis does not satisfy the 
Anderson condition. 

The band structure and DOS of the COO-IV state for 
the Pmca symmetry are shown in Fig. 11. For the Pmca 
symmetry, the band gap at the Fermi energy does not 
open; it is closed at the Y(0, i, 0) point and mostly closed 
at the X(i,0, 0) point. The band gaps at the L points 
which are on the direction from the centre of the FC4O4 
cube to the electron-rich sites are small in the COO-II 
state, and corresponding bands at the L(|, i, i) point of 
the two COO-II micro-layers come to the X(i, 0, 0) point 
and those of the L(— i, i, i) and L(i, — i, i) points are 
positioned at the Y(0, |,0) point. In contrast to the re- 
sults of the calculations with the Pmca lattice symmetry, 
the COO-IV state obtained assuming the Pmc2\ symme- 
try has the band gap at the Fermi level, and its magni- 
tude is - 2.5 x 10~ 2 eV for V = 0.3 eV and t pd = 0.35 eV. 
However, the band gap does not open for the small values 
of V. 

A V2a c x \/2a c x 2a c supercell comprising 64 B sites, 
which is the same dimension to the Cc symmetry cell, 
is also examined using the hopping integrals obtained 
from the Pmc2\ atomic coordinates. However, the COO- 
IV state is obtained as the ground state, and no further 
COO state is found. In this supercell, the X and Y points 
of the Pmca or Pmc2\ cells correspond to the T point, 
and the energy gap at the Fermi level is most narrow at 
this symmetry point in the first Brillouin zone. 

6. Discussions 

The electronic state in the high-temperature phase 
is still in controversy. Anderson proposed the state in 
which two Fe 3+ and two Fe 2+ ions are equally contained 
in every B site tetrahedron with the short-range order- 




zr x a y r do 12 



DOS (states/ eV atom) 

Fig. 11. Band structure and DOS of the COO-IV state obtained 
from calculations with V = 0.3 eV, i p( j = 0.35 eV and the hop- 
ping integrals evaluated from the Pmca atomic positions. The 
symmetry points Z(0,0,|), X(|,0,0), A(|,|,0), Y(0,|,0) and 
D(1,0, -i) with units 12L = ^2L = av^vr and 2jL = JL arc 
chosen. The dotted line represent the Fermi level. The band gap 
at the Fermi level does not open. 

ing. Yamada et al. have suggested a molecular polaron 
model. 50,51 On the other hand, half-metallic state has 
been predicted in LSDA calculations, 28,52,53 and results 
of resonant x-ray scattering measurements support the 
itinerant-electron picture. 54 

From the results of the HF calculation with the cu- 
bic structure, the COO-I state is anticipated to realise in 
the high-temperature phase. This is consistent with the 
fact that the magnetic easy axis for the high-temperature 
phase is in the (lll) c direction. However, the rhombo- 
hedral distortion expected in the presence of the COO-I 
state is not observed. Yamada et al. proposed a molecular 
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polaron model to explain the neutron diffuse scattering 
found in the high-temperature phase. 50,51 In this model, 
the valence fluctuation of F84O4 cubic cluster is strongly 
coupled with dynamical lattice distortion, i.e., the T<i g 
mode, within the cluster. However, the local lattice dis- 
tortions assumed in the model (see Fig. 1 in ref. 50) are 
the same as those in the presence of the COO-I and COO- 
II states as we discussed in § 5.2. Thus the diffuse scatter- 
ing can be interpreted as a result of the fluctuation of the 
fictitious orbital moment on the Fe404 clusters strongly 
coupled with dynamical lattice distortion. For this rea- 
son, wc consider the COO-I state remains short-range 
order in the high-temperature phase and thus no statical 
rhombohedral distortion occurs in the high-temperature 
phase. 

As for the low-temperature phase, variety of charge 
ordering model have been proposed. 9 ' 55 ' 56 However, the 
NMR 57 and x-ray anomalous scattering 58 cast doubts 
on the existence of charge ordering. Later on, mod- 
erate charge orderings with apparent local charges of 
Fe 2 4+ and Fe 2 6+ are inferred from the x-ray and neutron 
diffraction experiments. 3,4 The COO-IV state obtained 
for the low-temperature phase is also accompanied by 
charge disproportionation of Fe 26+ and Fe 2 4+ . However, 
the pattern of charge ordering predicted from the present 
theory is not the same to those inferred from x-ray and 
neutron diffraction experiments. 

We have described the COO-IV state as an alterna- 
tive stacking of microscopic layers of the COO-II state 
with two different orientations. However, the state can 
also be regarded as a COO-I state with the [lll] c mag- 
netisation direction modulated by the lattice distortion: 
the orbital magnetisation direction is slanted from direc- 
tion [lll]c to [011] c in the layer 1/4 < z m < 3/4 and 
direction [lll] c to [101] c in the layer -1/4 < z m < 1/4. 
Indeed, the COO-IV state can even continuously trans- 
form into the COO-I state as the values of t p d is increased 
from t p d = 0.5 eV with V = eV as demonstrated in 
Fig. 12. As was discussed in § 5.2, the COO-I state is 
expected to be accompanied by the rhombohedral cell 
distortion. This explains why the distortion of the unit 
cell in the low-temperature phase can be approximated 
by the rhombohedral distortion. 

As was mentioned before, the Jahn- Teller like distor- 
tions of oxygen octahedra surrounding the Fe ions at 
the Bl sites are observed in the experiments. 4 This dis- 
tortion found in the low-temperature phase can be ac- 
count for the presence of the COO-IV ordering as fol- 
lows. The Fe ions at the Bla and Bib sites and the 
neighbouring oxygens on the z m = 1/2 plane are il- 
lustrated on the left side of Fig. 13. The Fe-0 bonds 
along the cc-axis are elongated, and those along the y- 
axis are contracted at the Bla site. The Fe-0 lengths 
along the two directions are reversed at the Bib site. 
The complex orbitals at the Bla and Bib sites are \x~) 
and | Oil) in the COO-IV state (see Table V) and their 
electron clouds are stretched along the x- and y-axes, re- 
spectively. This antiferro-orbital alignment along the Bl 
chain is, therefore, stabilised by the Jahn- Teller like dis- 
tortion with alternative directions of the elongated Fe-0 
bonds along the chain. The same holds for the Bl chains 
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Fig. 12. Continuous transformation of the COO-IV state (a) into 
the COO-I state (b) with increasing value of t pc i from t pc i = 
0.5 eV. Alignments of orbital moments within —1/8 < z m < 1/4 
are shown. 
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Fig. 13. Illustration of the Jahn- Teller-like distortions along the 
Bl chain at z m = 1/2 (on the left) and the B3/B4 chain at 
2m = 5/8 (on the right). The pseudo- mirror plane perpendicular 
to the a m -axis presumed in the crystal structure analysis is also 
shown. The open and closed circles represent oxygen and Fe ions, 
respectively. The arrows from Fe ions represent the elongated Fe- 
O bonds. The Jahn- Teller-like distortions along the Bl chain are 
observed in the experiments, 4 and those along the B3/BA chain, 
which break the mirror symmetry, are also anticipated in the 
presence of the COO-IV state. 



on the z m = plane with the |101)/|y _ ) antiferro-orbital 
alignment. 

The similar antiferro-orbital alignment of |011)/|a; + ) 
or |y + )/|101) are present along the B3/B4 alternative 
chains parallel to the 6 m -axis in the COO-IV state. 
The alternative Jahn- Teller-like distortion is anticipated 
along the B3/BA chain from the theory. However, such 
distortions are still not found in the x-ray or neutron 
diffraction experiments. This chains are on the pseudo- 
mirror plane perpendicular to the a m -axis presumed in 
the crystal structure analysis (see Fig. 13) and search for 
such distortion on the B3/BA chains can be a check for 
the validity of the present theory. 

From results of the high-resolution photoemission 
spectroscopy (PES) experiments, 59 the insulating gap 
with the magnitudes ~ 0.1 eV and ~ 0.15 eV are inferred 
above and below Ty, respectively. Other PES experi- 
ments 60 and oxygen K-cdgc absorption measurements 61 
also support non-closure of the gap above Ty and in- 
dicate small-polaronic behaviour of charge carriers. The 
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formation of the pseudo-gap with negligible Drude weight 
in Ty < T < 300 K and a clear opening of the optical gap 
of ~ 0.14 cV below Ty are observed in measurements of 
the optical-conductivity spectra. 62 The COO states ob- 
tained with the cubic structure are either a semimetal 
with extremely low carrier density or a zero-gap semi- 
conductor and have a pseudo-gap with the magnitude 
~ 0.2 eV in the vicinity of Fermi level. In addition, the 
COO state is expected to be short-range order due to 
the strong coupling with the dynamical lattice distortion 
in the high-temperature phase. Thus, small-polaronic or 
incoherent-metallic behaviours of the high-temperature 
phase inferred from the above experiments are consistent 
with the present theory. 

Since the energy gap at the Fermi level is almost 
opened in these COO states obtained with the cubic 
structure, symmetry lowering caused by the monoclinic 
lattice distortion in the low-temperature phase would re- 
sult in the formation of the insulating gap. Indeed, the 
COO-IV state has a insulating gap of <~ 0.02 eV with the 
Pmc2i symmetry and large value of V ~ 0.3 eV. Inter- 
estingly, the absence of the inversion symmetry leads to 
the COO-IV state being multiferroic. 45 For the Cc super- 
cell, the band gap at the Fermi level is narrow or closed 
only at the T point. The size of this gap might be larger, 
if the rhombohedral cell distortion in combination with 
the Cc atomic displacements are further considered. 

7. Conclusions 

In this paper, a spinless three-band Hubbard model for 
tig electrons with minority spin on the B sites is inves- 
tigated with the Hartree-Fock (HF) approximation and 
exact diagonalisation method to discuss the electronic 
structure and Verwey transition in magnetite. 

Complex-orbital ordered (COO) states with non- 
collinear orbital moments ~ 0.4 are found as the 
ground state with the cubic structure of the high- 
temperature phase when the electron hopping via ad- 
jacent O 2p orbitals is strong. Two kinds of COO states 
with total orbital moment along directions (Hl) c (COO- 
I state) and (110) c (COO-II state) are obtained. These 
states are semimetal or zero-gap semiconductor, and 
pseudo-gaps of ~ 0.2 eV are formed near the Fermi level. 

The formation mechanism of the COO state is dis- 
cussed analytically and the stability of the state against 
the electron correlation effects is examined with the exact 
diagonalisation method for the 16-site system. The COO 
state does not stem from the orbital moments induced 
by ordered spin moments via the spin-orbit interaction 
as usual magnetic materials. It is an exotic magnetic or- 
dering within the orbital degree of freedom. The system 
of electrons with minority spin on the B sites can be 
described as a fee lattice of interacting Fe 4 tetrahedral 
molecules. From this viewpoint, the COO state originates 
from the formation of Hund's second rule state in spin- 
less pseudo-d symmetry molecular orbitals on each of the 
FC4 tetrahedron and ferromagnetic alignment of their fic- 
titious orbital moments. 

For the low-temperature phase, the HF calculations 
were performed with the hopping integrals evaluated 
from the Pmca and Pmc2i structural models. A COO 
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state, which can be described as an alternative stack- 
ing of microscopic layers of the COO-II states with two 
different (110) c directions of the magnetisation, is ob- 
tained (COO-IV state). The state is accompanied by 
moderate charge disproportion of Fc 2 6+ and Fe 2 4+ and 
in accordance with Jahn- Teller like distortion. Although 
we used the orthorhombic atomic positions, the symme- 
try of the COO-IV state spontaneously lowers to the 
monoclinic. Since the state can also be regarded as a 
modulated COO-I state with pseudo-trigonal symmetry, 
this explains conventional interpretation of the mono- 
clinic distortion: combination of the rhombohedral cell 
distortion and atomic displacements in the cell with the 
orthorhombic symmetry. 

From these findings, we interpret the Verwey transi- 
tion as follows. Above Ty, the COO-I state remains to 
be short-range order impeded by dynamical lattice dis- 
tortion. The strong coupling between the COO-II states 
and lattice leads to the formation of the COO-IV state 
below Ty. 
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Appendix: MOs in Isolated Tetrahedron 

The four ti g orbitals of the isolated tetrahedral 
molecule can be reduced to five different basis sets of the 
irreducible representations, i.e., A\ +E + T\ + 2T 2 in the 
Td point group symmetry. The energies s(T) and bases 
|r, 7) of these irreducible representations T are listed in 
Table A-l. There are two different basis sets for the T2 ir- 
reducible representation and here, we call them T£ and 
T£. The hopping integral between them (T?\H\T£) is 
also shown. For each of the basis sets of the three dimen- 
sional representations T\ and T2, only one of the bases is 
presented; the other two bases, i.e., those with /? and 7 
for T\ and 77 and ( for T2, are easily obtained from the 
bases in the table by rotating ±27r/3 around the [lll] c 
axis of the tetrahedron. 
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Table A-l. Energies e(T) and bases |T,7) of irreducible represen- 
tations r in the symmetry on the t^g orbitals of tetrahedral 
molecule. 
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